Universal properties of Fermi gases in arbitrary dimensions 
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We consider spin- 1/2 Fermi gases in arbitrary, integer or non-integer spatial dimensions, interact- 
ing via a Dirac delta potential. We first generalize the method of Tan's distributions and implement 
short-range boundary conditions to arbitrary dimension and we obtain a set of universal relations 
for the Fermi gas. Three-dimensional scattering under very general conditions of transversal confine- 
ment is described by an effectively reduced-dimensional scattering length, which we show depends on 
the three-dimensional scattering length in a universal way. Our formula for non-integer dimensions 
interpolates between the known results in integer dimensions 1, 2 and 3. Without any need to solve 
the associated multichannel scattering problem, we find that confinement-induced resonances occur 
in all dimensions different from D = 2, while reduced-dimensional contacts, related to the tails of 
the momentum distributions, are connected to the three-dimensional contact by a correction factor 
of purely geometric origin. 
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I. INTRODUCTION 

The dependence of physical properties on dimensional- 
ity is an important issue with great theoretical and prac- 
tical ramifications. The most well-known examples are 
arguably theories containing extra dimensions, such as 
the original attempt of Nordstrom, Kaluza, and Klein 
to unify gravity and electrodynamics by adding an ex- 
tra spatial dimension |l|. These ideas are currently a 
backbone in modern thereotical high-energy frameworks 
However, the quest to go in the other direction, 
i.e. to systems with spatial dimension less than three, is 
equally important since low-dimensional structures such 
as wires and surfaces dominate our daily surroundings. 
In fact, quantum dynamics in one- and two-dimensional 
systems can be very different from the three-dimensional 
case 0. 

An obvious mathematical question related to variation 
in dimensionality concerns interpolation between integer 
values, and its relevance in physics. This notion, however, 
turns out to have great value in such areas as the charac- 
terization of chaotic phenomena Qjin surface growth 8] , 
on the study of phase transitions [9| , for domain forma- 
tion in spin models [l^l, and also for quantum computa- 
tion [ll| on fractal lattices [l^ • From a more formal the- 
oretical point of view, variation of dimensionality beyond 
integer has had tremendous impact on particle physics, 
statistical physics, and condensed-matter through the fa- 
mous e-expansion by Wilson [isj and dimensional regu- 
larization by t'Hooft [l^. 

The advantage of studying cold atomic gases is the 
clean conditions that allow an isolation of specific effects 
of interest 15]. The experimental control over the inter- 
actions means that one has direct access to the strongly- 
coupled regime where the scattering amplitude in a two- 
component Fermi gas saturates the unitarity bound of 
quantum mechanics (or equivalently when the s-wave 
scattering length diverges). An experimental surge has 



managed to map out important features of the univer- 
sal thermodynamics of this system [l^. However, as 
shown by Tan [Tt'I (see also [Ts'l), the unitary regime 
has a direct connection between microscopic quantities 
such as the momentum distribution and the macroscopic 
thermodynamic observables. These relations have sub- 
sequently been demonstrated in experiments on three- 
dimensional gases [l^ll^l- In lower dimensions, similar 
microscopic relations have been derived [2T| - [26| . and ex- 
periments testing these relations in two dimensions are 
on-going 

In the current paper we provide an interpola- 
tion/extrapolation formula for Tan's contact relations as 
a function of the spatial dimension, D, which is consid- 
ered a continuous variable in the general setup that we 
discuss. The objective is to provide interpolation for- 
mulas for the contact identities that can be used for in- 
stance in interpolating from known limits into unknown 
territory, and testing the results of different analytical 
and numerical approaches to strongly-coupled dynam- 
ics. This is similar in spirit to the e-expansion based 
on the properties of the Fermi gas in = 2 and D = A 
interpolated to £> = 3 and compared to Monte Carlo 
calculations by Nishida and Son [28|, [2^. These rela- 
tions should be useful in general strong-coupling investi- 
gations, for instance when applying the renormalization 
group using the e-expansion. More generally, they may 
provide hints of whether certain confining geometries al- 
low access to systems where the dynamics indicate an 
effective non-integer dimensionality. An example could 
be a holographic model where the boundary turns out to 
have fractal dimension, as for instance seen in the spin 
models mentioned above. 

Using these universal relations wc find that, under 
very general conditions, effective interactions in quasi- £)- 
dimensional Fermi gases, which appear as emergent theo- 
ries for transversely confined three-dimensional systems, 
are related to the original three-dimensional interaction 
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through a universal functional form. This new type of 
emergent universality has a strong predictive power: the 
so-called confinement-induced-resonances [sOl are shown 
to occur in all dimensions except for the two-dimensional 
case, which a gree s with the known results in integer di- 
mensions [30II3H. 

The paper is structured in the following way. In Sec. 
im we introduce Tan's ?7-distribution in terms of inte- 
grals over D-dimensional momentum space. In Sec. IIIIl 
we present the corresponding set of Tan relations, pro- 
viding connections between a number of microscopic and 
macroscopic properties of the Fermi gas in D dimensions. 
In Sec. IIVI we recall Olshanii's and Petrov-Shlyapnikov's 
results for scattering confined to one and two dimensions, 
respectively, and we show that they are special results 
of Z?-dimensional expressions, which link the scattering 
lengths and contacts to their three-dimensional equiva- 
lents. Sec. fVl concludes the manuscript. 



II. FORMALISM 

We consider a spin- 1/2 Fermi gas with contact s-wave 
two-body interactions in arbitrary spatial dimension D. 
The scattering of a pair of fermions is then completely 
characterized by a single quantity, namely the scattering 
length Contact interactions have to be regularized in 
D > 2 and renormalized (in any dimension) to fit the cor- 
responding D-dimensional scattering amplitude for two- 
particles or equivalently, when a bound state exists, to 
the two-body binding energy in vacuum. This proce- 
dure can be done in a number of ways, all equivalent 
to each other2_ non-perturbative Gell-Mann-Low renor- 
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malization [32l . l33l |. dimensional regularization 
i'-dimensional Bethe-Peierls boundary conditions 
Fermi-Huang pseudopotential [37i | or its A- fam ily |38| and 
the more recent method of Tan's distributions [l7ll22ll39l - 
IZlj . In this work, we will make use of the latter method. 

In order to treat non-integer-dimensional spaces and, 
in particular, integrals of functions of £)-dimensional 
variables, we first need to define Euclidean space in ar- 
bitrary dimension. There are several ways of doing this, 
e.g. Stillinger [i^ and Wilson spaces Q. The spaces 
differ in the fact that Stillinger spaces do not have vec- 
torial character, while Wilson spaces do. Integration in 
these two spaces can be successfully defined, although 
the appearance of negative integration weights seems un- 
avoidable. For all purposes of this paper, both Stillinger 
and Wilson spaces are equivalent, and we will only need 
simple integrals where no conceptual problems arise. In 
particular, for (hyper-)spherically symmetric integrands 
we can safely use the replacement: 
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where \{D) = 27r^/2/r(D/2) is the D-dimensional sohd 
angle. 



We now derive Tan's rj distribution 17[ in arbitrary 
non-even dimension D. Note that the case D = 4 is not 
relevant — there is no scattering by a Dirac delta even 
after renormalization — while for D — 2 we derived the 
corresponding expression in [2^. The rj distribution is 
defined via ri{k) = 1 (fc < 00) and 
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where is given by [44 



r2(i:)/2)sin(7r(D-2)/2) 



tt{D-2)/2 
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(2) 



(3) 



and coincides with the two-body binding energy if the 
scattering length ao is positive. The ?7-distribution, as 
defined in Eq. ([2]), provides a simple way of implement- 
ing short-range boundary conditions for a D-dimensional 
delta interaction and its use is completely equivalent 
to cut-off and dimensional regularization-renormalization 
procedures |4l|. 

From previous analyses for dimensions D = 1,2,3 [H, 
m. Ell, we see that 



r^{k) = l-g{k)5{l/k), 



(4) 



where Q{k) is determined from Eq. ([2]) and is therefore 
given by 
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The series representation of the above integral for non- 
even dimensions is obtained from 
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We note that terms higher than n — [-D/2J in the series 
([7]) do not contribute to the integral in Eq. ([2]). 



III. TAN RELATIONS IN D DIMENSIONS 

We show here the set of universal relations for Fermi 
gases in arbitrary dimensions. We single out the two- 
dimensional case, for obvious reasons, which has been 
derived previously using Tan's distributions [22l |. Pro- 
ceeding in a way which parallels refs. [U H^], we obtain 
the complete set of Tan relations. 
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a. Energy theorem The energy Eo satisfies the fol- 
lowing energy theorem 



(8) 



where W is a single particle trapping potential acting 
on each fermion in the gas and nk,o- is the occupation 
number of the momentum state with momentum k and 
spin a. Using Eqs. (I4I5I6I7[) . we can expand the above 
relation to obtain 
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where Vd is the _D-dimensional hypervolume and Cd — 
linik-j-oo If^nk.CT is called the ZJ-dimcnsional contact. 
b. Adiabatic relation We obtain 



dEp _ apCp 
dap On"^ 



(10) 



where 



ap = {D- 2) 



2\ 2--D/2 



VDXiD)Tr 



csc{^D/2){a]jEB)°'^-^. 



(11) 



Generalized virial theorem Given a trapping po- 



tential of the form W (x r^, we obtain 
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d. Pressure relation If the system is homogeneous, 
that is, W = 0, the pressure relation reads 



PnVp = 1 



2Ep 



apCp 

„D-2 



(13) 



e. Two-particle loss rate Adding a small imaginary 
part to the scattering length, that is, replacing ap with 
ap + iai, we obtain the last of the relations for the two- 
body loss rate Tp 
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IV. 



EFFECTIVE D-DIMENSIONAL 
INTERACTIONS 



The effective reduction of a three-dimensional zero- 
range interacting two-body system placed under the 



influence of a (3 — D)-dimensional trap to an inte- 
ger lower dimension {D — 1,2) is achieved usin g a 
method due to Olshanii ^30] or variations thereof [3lll45l. 
l46j . For the effective one-dimensional case under two- 
dimensional isotropic harmonic confinement the reduced 
one-dimensional scattering length is given by f30l | 



ai = h 

03 V2 



(15) 



where £q = ^/hfrrujj is the oscillator length, x = 
Vmis^ooU^ dtt-^/^ - E^'=is'"'/') ~ 1-46. In two di- 
mensions, under a one-dimensional transversal harmonic 
trap, the effective scattering length reads [3lj 



a2 = 2e-^^W-£oexp(-vW24/a3), 



(16) 



where is Euler's number and _B « 1. 

Actually, dimensional reduction can be done from any 
D' 7^ 3 to a lower dimension D, but we will assume that 
the reference space has dimension three. The general idea 
also applies to non-integer D and involves the calculation 
of zero-energy (up to a constant offset) wave functions ^ 
that behave as 



*(r) ^ V(rD)0o(r3-D), 



(17) 



as rp — > CX3, where vp represents the Z)-dimensional rel- 
ative position, ip (rp) is the _D-dimensional zero-energy 
wave function corresponding to an effective scattering 
length ap (to be calculated), and represents the 
ground-state mode of the (3 — Z?)-dimensional external 
trap. 

This kind of effective dimensional reduction is a mul- 
tichannel scattering problem with an infinite number of 
coupled channels and, obviously, requires regularization 
and renormalization due to the three-dimensional zero- 
range interaction. The multichannel problem needs to 
be solved, in principle, in order to obtain the functional 
form of the I?-dimensional scattering length a^i in terms 
of the original three-dimensional one 03 in vacuum, which 
can be regarded as the renormalization condition. 

We here derive the general functional form for ap = 
F{a3,D) using only dimensional analysis and the adi- 
abatic theorem, Eq. (fTO|) . as follows. We first need a 
length scale, which we denote by xp which, in the case of 
an isotropic harmonic confinement must be proportional 
to the oscillator length. The I?-dimensional contact and 
energies are denoted by Cp and Ep, respectively. The 
number of identical reduced-dimensional subsystems em- 
beded in the three-dimensional space is denoted by Afp 
(typically this is realized through confinement in optical 
lattice potentials such that in D = 1 each of these 
systems is a tube while in D = 2 they are pancakes) . As- 
suming that each of these subsystems is well separated 
from the rest and has energy Ep, the total energy of 
the system E3 will be given by E3 = NpEp, up to an 
as-independent energy offset. The three-dimensional vol- 
ume of each subsystem is given by its hypervolume Vp 
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times its transversal extent ^, so the total volume 
of the system is given by V3 = Afox^^Vn- Using these 
facts in the three- and D-dimensional adiabatic relations, 
Eq. ([TU|) . we obtain 



oidCl 



dF 



D 



D-l 



das 



{a3,D) 
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(18) 
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From Eqs. (US]) and (HH), we see that F'^^^dF/d, 
a'^^ in one and two dimensions. Assuming that this is 
also the case in arbitrary dimension D, we arrive at the 
condition 



dF 



[F{a3,D)]D-ida. 



(19) 



where no is independent of 03. From Eq. (jl9l) and for 
D ^ 2 we obtain 



an = F[a3,D] = 



l/(D-2) 



(20) 



_{D-2)i-fD + Ko/aa) 
where 7c is a constant. In two dimensions we obtain 

a2 = F[a3, 2] = e'^^e"''^/''^ (21 



The above expressions ([201) ^-^d (|2T|) are the desired rela- 
tions between effective Z?-dimensional scattering lengths 
and the three-dimensional scattering length in vacuum. 
For all dimensions D there are two a priori unknown con- 
stants and Kd. To obtain the constant we need 
to solve the full multichannel scattering problem as dis- 
cussed above, but ko can be exactly calculated from a 
mean- field shift in the coupling constant, as shown in the 
Appendix. Note that, as must obviously be, Eq. (f20|) for 
D = 1 and Eq. (f2T|) agree with the known results (fT5|) 
and (fTB]) . respectively. The case D = 3 is trivially satis- 
fied, with K3 = 1 and 73 = 0. 

The above analysis has some striking, non-trivial con- 
sequences. First, one obtains a confinement-induced res- 
onance (CIR) corresponding to a^i = for all < 2 
whenever 03 = —kd/id- For all dimensions 2 < D < 3 
this matching condition results in od — >■ 00, and it is 
only the two-dimensional case where no resonant behav- 
ior may occur. Second, while the values of k^i and are 
mo del- dependent, relations ((20l) and (|211) are universal. 
These are valid for any isotropic transversal (3 — D)- 
dimensional confinement potential provided that either 
(i) (3 — £')-dimensional center-of-mass and relative coor- 
dinates can be separated [13] (as is the case for harmonic 
confinement) or (ii) the corresponding CIR is so broad 
as compared to further resonances that a so-called single 
pole approximation [i^ for the effective scattering length 
is valid, which is typically the case [i^ . 

If we use now Eqs. (HI]) and (HH) in Eq. ([H]), in the 
case of isotropic harmonic confinement we find the pro- 
portionality constants that relate the effective contacts 



in D = 1, 2, to the three-dimensional one: 

1 X2 



C2 = 



27r 4 



C3 



(22) 
(23) 



For general dimension, we can obtain the above relation 
up to a multiplicative constant and it reads 



C3 



(24) 



where S^jj — ajj/Vo- Since we have cx ft^/m, we have 



from dimensional grounds that kjj 



n3-D 



and Ac is 



dimensionless and a3-independent. Therefore, A^i is of 
purely geometric origin and we call it the D-dimensional 
geometric correction factor (GCF). 

Some important remarks about the GCF are in or- 
der. We have assumed that the system is tightly con- 
fined by a (3 — D) dimensional trap, which implies that 
the effective D-dimensional Fermi energy is much smaller 
than the first excited energy (oc fiw) of the transversal 
trap, E^j^^ <^ hw. The GCF is in this case a constant, 
but the D-dimensional and three-dimensional contacts 
are not necessarily the same. To sec this, we can esti- 
mate the two- and one-dimensional transversal extents 
X2 and x\ , respectively, for the case of harmonic confine- 
ment. We define (r|,) — {tpo\r'^ where jf/'o) repre- 
sents the single-particle ground-state in the transversal 
trap. Then, the quasi-two-dimensional system will have 
a extent given by X2 — 2^ (rj) = v^^Oi which gives 
A2 = I/Vtt- In the one-dimensional case, assuming a 
cylindrical shape for the quasi-one-dimensional tubes, we 
have a;^ = 7r(r|) = tt^q, which gives Ai = 1/2. These cor- 
rection factors must be taken into account when relating 
C3 to the Z?-dimensional contact Cd, in order to ver- 
ify universality in effectively reduced-dimensional Fermi 
gases. 

In a recent experiment concerning CIRs (49| . anhar- 
monicities in the trapping potential are present, which 
are not described by our previous analysis. This im- 
plies that close to anharmonic CIRs [l^jSgl, the quasi- Z?- 
dimensional contact Cd is not proportional to the three- 
dimensional contact C3. In the quasi-lD isotropic case 
we may, on the other hand, use a two-pole approach - 
see ref. (45j for the single-pole approximation - to de- 
scribe the presence of a second resonance, which can be 
justified in the spirit of Fade approximants. We assume 
that the quasi-lD interaction strength gi = -~2h^ /mai 
has the form 



1 



1 



51 



ri + rjia3 r2 + r22a3' 



(25) 



where F^ and fij (i = 1,2) mark the positions Og — 
—Ti/Vti of the two CIRs, while Vti/Q,2 defines the relative 
width of the resonances when they are well separated 
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,(2) I 



^1 for i = 1, 2). We then have, close 



to the broad CIR (e.g. 03 



C; 



Ci ill 



^3 ^ 1 I OT^ ^ I ^1 



from Eq. (fT9)) . 



- 1 



(26) 



The above expression gives the hnear correction to the re- 
lation between 3D and quasi-lD contacts in the isotropic, 
anharmonic case. The first term on the r.h.s. of Eq. ((26)) 
corresponds to the proportionality constant in the har- 
monic case, while the linear term reflects slight deviations 
due to anharmonicity in the vicinity of the broad CIR. 



onto lower dimensional physical realization of the pre- 
dicted dynamics. These theories are used in addressing 
strongly-coupled dynamics, which is a limit in which Tan 
relations have proven particularly powerful in recent cold 
atomic gas experiments. We hope that the current work 
can provide a foundation for exploring the implications 
of Tan relations within other fields as well. 
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V. CONCLUSIONS 

We have derived universal Tan relations for spin-1/2 
Fermi gases interacting via a Dirac delta potential in ar- 
bitrary integer and non-integer spatial dimensions. The 
formalism introduced was examplified by considering ex- 
trapolation of known results for the three-dimensional 
case where the contact of the unitary Fermi gas has 
been calculated for the ground-state energy using the e- 
expansions ^2^ [2^ . Starting from this three-dimensional 
initial condition, we studied the properties of Fermi 
gases subject to transversal (3 — £')-dimensional con- 
finement which results in an effective realization of D- 
dimensional Fermi gases. Using the universal relations, 
we found that, in the limit of tight transversal confine- 
ment, the effective D-dimensional two-body interactions 
exhibit a universal functional dependence on the three- 
dimensional scattering length in vacuum. This universal 
function implies the appearance of confinement-induced 
resonances [30l . [sij in any integer or non-integer dimen- 
sion, with the only exception being the well-known two- 
dimensional case. Furthermore, the quasi- D-dimensional 
contact is related to the three-dimensional one, by a cor- 
rection factor whose value is dictated by the geometry 
of the transversal trap. The example can be general- 
ized to the case of a _Do-dimensional reference space and 
{Dq — Z?)-dimensional confinement, such as is often seen 
in higher-dimensional theories {Dq > 3) where our every- 
day observable Universe is embedded as a submanifold. 
The prescription provided here shows how to extrapo- 
late higher-dimensional universal quantities down to the 
relevant submanifold. 

Tan relations are of relevance for the e-expansion (l3| 
around two and four dimensions, and can be of use in 
other setups where dimensional interpolation or extrap- 
olation of physical quantities is necessary. A case of in- 
terest could be extradimensional or holographic theories 
where higher-dimensional observables must be projected 



Appendix A: Calculation of kd through the 
mean-field shift in the coupling constant 

In general, the effective Z3-dimensional coupling con- 
stant is related to the scattering length via 



go = {D - 2)\{D){n^ /m)a\ 



(Al) 



When the system is transversely confined by a (3 — D) 
dimensional trap, it is given by [30| 



=53|0o(r3-D = 0)|2 + O(52), 



(A2) 



to lowest order in 173, where <f>^ is the normalized ground- 
state wave function of the (3 — D)-dimensional trapped 
system in the relative coordinate. 

For the case of a harmonic trap, we obtain (soj 



io(r3-Z5 = 0)p 



1 



(2^)(3-D)/2^3-C- 

Equation (|A2[) becomes in this case 

53 



9d = 



(27r^§)(3-^)/2 



(A3) 



(A4) 



Using the corresponding expression for gjj in terms of 
a£), we find 



An 



{D-2)X{D){2n£l)(^~D)/: 



l/(D-2) 



(A5) 



to lowest order. From Eq. pO|) to lowest order in 03 we 
then obtain 



Kd 



2(D-i)/2Y(D/2) ° 



(A6) 
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